Temperature dependence of volume and surface symmetry energy coefficients of nuclei 
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The thermal evolution of the energies and free energies of a set of spherical and near-spherical 
nuclei spanning the whole periodic table are calculated in the subtracted finite-temperature Thomas- 
Fermi framework with the zero-range Skyrme-type KDEO and the finite-range modified Seyler- 
Blanchard interaction. The calculated energies are subjected to a global fit in the spirit of the 
liquid-drop model. The extracted parameters in this model reflect the temperature dependence 
of the volume symmetry and surface symmetry coefficients of finite nuclei, in addition to that of 
the volume and surface energy coefficients. The temperature dependence of the surface symmetry 
energy is found to be very substantial whereas that of the volume symmetry energy turns out to be 
comparatively mild. 

PACS numbers: 21.10.Dr, 21.30.Fe, 21.65.Ef, 26.30.-k 
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The symmetry energy coefficient a v sym of infinite nu- 
clear matter is conventionally defined by the relation 
e(X) = e(X = 0) + a v syrn X 2 . Here e is the energy 
per nucleon of the system at isospin asymmetry X = 
(p n — Pp)/(p n + pp), p n and pp being the neutron and 
proton densities, respectively of the system. For homo- 
geneous nuclear matter, this definition works extremely 
well, e(X) is seen to be bilinear in X for nearly all 
values of asymmetry [l], For warm nuclear matter, 
the symmetry free energy coefficient fg ym is likewise ob- 
tained from f{X, T) - f(X = 0, T) = f" ym {T)X 2 , where 
f(X,T) is the per-nucleon free energy of the matter at 
asymmetry X and temperature T. These asymmetry co- 
efficients are measures of the energy or free energy release 
in converting asymmetric nuclear system to a symmetric 
one. For infinite nuclear systems at saturation density 
po and temperature T = 0, the value of a v sym is usually 
taken in the range of - 30-34 McV 041. 

In the global fitting of the nuclear masses in the frame- 
work of the liquid-drop mass formula, the symmetry co- 



efficient 



< sym 



enters as a phcnomenological parameter. 



Nuclei being finite systems, it is realized that varying 
density profiles of different nuclei necessitate introduc- 
tion of a mass-dependent surface component in a sym (A) 
in addition to the mass-independent volume component 
a v sym . In the literature, two different definitions have 
been used for a sym (A). The first, hereafter referred to as 



is. 
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1 + 

PE 

and the second, hereafter referred to as II @ is, 
Qsym{A) = a sym — a sy m A ^ ■ 



(1) 



(2) 



and that of a v sym / '(3e is in the close range 
2.4±0.4 



In definition I, /3g is a measure of the surface symmetry 
energy, a s sym is the surface symmetry energy coefficient in 
definition II. In the limit of very large A, {a v sym ) 2 / '/3 'e ~ 

is taken as 

45 MeV " 

Of' 

It is evident that the symmetry energy coefficient has 
an extremely important role in describing properly the 
nuclear binding energies along the periodic table and in 
getting a broad understanding of the nuclear drip lines. 
It also plays a seminal role in guiding the dynamical evo- 
lution of the core collapse of a massive star and the asso- 
ciated explosive nucleosynthesis. A large (small) magni- 
tude of a sym inhibits (accelerates) change of protons to 
neutrons through electron capture El • This change 
in isospin asymmetry has its import in the nuclear equa- 
tion of state (EOS) and thus on the dynamics of the 
collapse and explosive phase of a massive star. Matter in 
that phase is warm, it is therefore essential to know with 
precision the thermal dependence of the symmetry coeffi- 
cients. Furthermore, in this collapse or bounce phase, the 
nuclear matter is inhomogeneous; it nucleates to clusters 
of different sizes. Knowledge about the thermal evolution 
of the symmetry coefficients of finite nuclei then becomes 
a matter of central importance. 

In the low temperature domain (T <2 MeV), calcula- 
tions of the symmetry coefficients of atomic nuclei have 
been done earlier by Donati et.al [l2| in a schematic 
model. The motion of the nucleons in a fluctuating mean- 
field results in a nucleon effective mass that carry signa- 
tures of nonlocality in space (the fc-mass m^) and also 
nonlocality in time (the energy-mass m u ). The en ergy 
mass m w is seen to decrease with temperature (T3 , Il4j . 
this brings in a decreased density of states and thus an 
increase in the symmetry coefficient. Calculations in 
this limited temperature range have further been done 



by Dean et.al [15[ in a shell model Monte-Carlo frame- 
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work. It provides qualitative support to these earlier find- 
ings. The symmetry coefficients, however, are found to be 
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of the temperature dependence of the volume and surface 
symmetry coefficients of nuclei have also recently been at- 
tempted by Lee and Mekjian [l6j in a density functional 
theoretic approach. These calculations are also limited 
to low temperatures (T <3 MeV); the approximations 
employed here keep the results meaningful in this small 
temperature domain. 

Exploring the thermal evolution of the symmetry co- 
efficients of specific atomic masses has been attempted 
fl7| in a broader temperature range (T < 8 MeV) more 
recently. The energies and free energies of the hot nuclei 
are calculated in the finite-temperature Thomas-Fermi 
framework (FTTF) with the subtraction technique [l8| 
with suitable choice of effective interactions. Dynamical 
changes in the energy-mass m w are taken care of. For a 
nucleus of mass A, the symmetry coefficient is defined as 



l syrn 



(A, T) = [e n (A, X u T) - e n (A, X 2 ,T)]/(X 2 - A 2 2 ).(3) 



Here e n 's are the nuclear part of the energy per nucleon 
of the nuclear pair of mass A but having different charges 
and X\ and X 2 are the asymmetry parameters of the nu- 
clei. For a finite nucleus with Z protons and ./V neutrons, 
X is defined as (N — Z)/A. Similar to a sym (A,T), the 
symmetry free energy coefficient f sym (A,T) can be de- 
fined. These definitions suffer from the fact that unique 
values of a sym or f sym for a nucleus of mass A can not be 
prescribed; the values depend on the choice of the isospin 
asymmetric nuclear pair. 

The present communication is aimed to arrive at un- 
ambiguous values of the temperature dependence of the 
symmetry coefficients. For a set (sixty nine) of spherical 
and non-spherical nuclei covering almost the entire pe- 
riodic table (we take 36 < A < 218 and 14 < Z < 92, 
the list of the nuclei is taken from Ref. [13); the en- 
ergies and free energies are calculated in the subtracted 
FTTF procedure, taking into account the dressing of the 
nucleon mass to energy-mass that arises from the 
coupling of the nucleonic motion with the surface vibra- 
tions [HI, [3 [2(| ■ Two effective interactions are chosen, 

i) the zero-range Skyrme-type interaction KDEO [2l[ and 

ii) the finite-range modified Seyler-Blanchard (SBM) in- 
teraction. The KDEO interaction reproduces the binding 
energies of many nuclei ranging from normal to exotic 
ones with a deviation which is much less than 0.5% for 
most cases. In addition, it has been extremely success- 
ful in reproducing the breathing mode energies of many 
nuclei, their charge radii and spin-orbit splitting. The 
SBM interaction also has been very successfully applied 
in getting properly the ground state binding energies [22| , 
charge rms radii, giant monopolc resonance energies etc. 
[H |24|. The SBM interaction is given by 



// ( r , P, P) = c l,u vi (r, p) + v 2 (r, p) 
«i = -(l-p 2 A> 2 )/(ri,r 2 ), 



TABLE I: The parameters of the KDEO effective interaction 

to tj t2 t-j, Xq 

(MeV-fm 3 ) (MeV-fm 5 ) (MeV-fm 5 ) (MeV-fm 3(1+a) ) 



.-i-i 



■f-2 



-2526.52 



430.94 



-398.38 



14235.52 



0.7583 -0.3087 -0.94 



TABLE II: The parameters of the SBM effective interaction 
(in MeV fm units) 



C, 



( u 



d 



348.5 



829.7 



0.6251 



927.5 



0.879 



1/6 
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/(ri,r 2 ) 
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The strength parameters Ci for like pairs (n-n,p-p) and 
C u for unlike pairs (n-p) carry information on the isospin 
dependence in the interaction. The densities at the sites 
ri and r 2 of the two interacting nucleons with momenta 
Pi and p 2 are given by p(vx) and p(r 2 ); r = |ri — r 2 | 
and p = |p x — p 2 |. The range of the interaction is a; b 
and d are measures of the momentum and density depen- 
dence in the interaction and k controls the stiffness on 
the nuclear EOS. The procedures for determining these 
parameters are given in detail in Refs. [U, [25|. The 
parameters for KDEO and SBM interaction are listed in 
Table I and II, respectively. The values of the satura- 
tion density p s , the volume energy, the isoscalar volume 
incomprcssibility A'oo, the volume symmetry coefficient 
a v sym , the symmetry incomprcssibility K sym , the symme- 
try pressure L and the critical temperature T c for these 
two interactions are listed in Table III. It is worthwhile 
to note that the values of the symmetry coefficients a v sym , 
K sym and L lie in the range suggested by the empirical 
constraints emerging out of the analyses of different re- 
cent experimental data [26T4291 ] . The method for obtain- 
ing the density profiles of hot nuclei and their binding 
energies in the subtracted FTTF approach, with subse- 
quent modification due to energy-mass with the SBM and 
Skyrme-type interaction has been described in some good 
detail in a recent article [lTj ; we therefore do not repeat 
it here. The energies and free energies of the chosen sixty 
nine nuclei arc calculated with this prescription in a tem- 
perature grid. At a particular temperature, the energies 
are then fitted in the framework of the Bethe-Weizacker 
mass formula 



E(N,Z,T) = a v (T)A + a s (T)A 



2/3 
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(A,T)X 2 A, 



(0) 



F(N7,T\ = f„.(T\A+fJT\A 2 / 3 + a. 



Z' 1 



3 



fsym(A, T)X 2 A. 



(7) 



In Eq. (6), and a sym are the volume, surface, 

Coulomb and symmetry energy coefficients. Similarly, 
f v ,fs and fsym are the corresponding free energy coef- 
ficients. The Coulomb energy has an implicit tempera- 
ture dependence; it does not contribute to entropy. Since 
they are precisely known in a calculation, we try to make 
a four-parameter fit with only the nuclear part of the 
energies and free energies, 

E n (N, Z, T) = a v (T)A + a s {T)A 2 ^ + a syrn (A, T)X 2 A,(8) 



F n (N, Z, T) = f v (T)A + f s (T)A 2 / 3 + f sym (A, T)X 2 A.(9) 

Here E n and F n are the nuclear part of the energy and 
free energy of the nucleus; a sym (A, T) is given by Eq. (I) 
or Eq. (2). In a similar spirit, f sym (A,T) is written as 



fsym (-^3 -^"0 



f sym (-^') 



/™m( T )^_l/3 



}sym{A,T) — fs ym (T) — fs ym {T)A 

The four-parameter set f v ,f s ,fg, 



-1/3 



(10) 

(11) 

(or p F ) 



'sym an( l fsyn 

have the same connotation as the set a v , a s , a v sym and 



a 



sym 



(or Pe), except that the former set refers to free 
energy. The parametric values of the volume energy a v 
and the volume symmetry free energy /„ are shown as a 
function of temperature in panels (a) and (b) of Fig. 1. 
At T=0, a v (or /„) very closely reproduces the energy 
per nucleon of symmetric nuclear matter. At low tem- 
peratures, a v and f v are nearly independent of the in- 
teractions chosen, at higher temperatures, a slight de- 
pendence is observed. For a particular interaction, these 
values, however, do not show any significant dependence 
on the chosen set I or II. Both a v and f v are seen to 
change quadratically with temperature. They are very 
well approximated with a v (T) = e(T = 0) + T 2 /A'i 
and f v (T) = f v {T = 0) - T 2 /K 2 , with K\ ~ 15.5 MeV 
and K<z ~ 24.0 MeV. It is to be noted that for infinite 
matter at a particular density and temperature, the en- 
ergy and free energy are canonically related (the entropy 
S = —(dF/dT) p , whence K\ = K 2 ); in the present case, 
density is a varying profile, also a v (T) and f v (T) are 
obtained from a least-squares fit to the energies of a mul- 
titude of nuclei. This may explain the different values of 
K\ and K 2 - 

In the right panels (c) and (d) of Fig. 1, the thermal 
evolution of surface energy and the surface free energy 
coefficients are shown. The surface energy (upper panel) 
increases slowly with temperature; with the KDE0 inter- 
action, a slight fall at very high temperatures is, how- 
ever, observed. With temperature, the surface free en- 
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FIG. 1: (color online) The volume energy and volume free 
energy coefficients shown as a function of temperature in the 
left panels (a) and (b). The full and dashed lines correspond 
to parametrization I and II, respectively. The blue color refers 
to the KDE0 interaction, the red color refers to SBM. In the 
right panels (c) and (d), the thermal dependence of surface 
energy and surface free energy is shown. 



TABLE III: Some bulk properties for infinite nuclear matter 
at the saturation density for the KDE0 and SBM effective 
interactions. The values of saturation density are in fm -3 
and all other quantities are in MeV. 



Force 


Ps 


a v 


K x 




a" 


L 


T c 


KDE0 


0.161 


-16.1 


229 


-144 


33 


45.2 


14.7 


SBM 


0.154 


-16.1 


238 


-101 


31 


59.8 


14.9 



different parametric forms for the dependence of surface 
free energy on temperature have been used. We find that 
the form of the type f s = f s (T = 0)[(T 2 ~T 2 )/(T 2 +T 2 )} a 
used in Ref. [31( gives a reasonably good fit with our cal- 
culated values for both the interactions using both the 
parameter set I and II with a ~ 0.95. 

In Fig. 2, the evolution of the volume symmetry energy 
a v sym and the volume symmetry free energy f v sym coeffi- 
cients with temperature are displayed in panels (a) and 
(b), respectively. The behavior of a v sym depends on how 



L sym. 



(A) is defined. In definition I, it falls with tempera- 
ture, in definition II, it shows a slow increase. The nature 
of the fall of a v sym (in I) or its increase (in II) is nearly 
the same for both the interactions. The coefficient fg ym , 
however, shows nearly no dependence on temperature for 
both the interactions and in both definitions. 

In the left panels of Fig. 3, the thermal dependence of 
the coefficients (3e and (3f as used in Eqs (1) and (10) 
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FIG. 2: (color online) The volume symmetry energy and the 
volume symmetry free energy coefficients plotted as a func- 
tion of temperature. The lines and the colors have the same 
meaning as in Fig. 1. 
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FIG. 4: (color online) The symmetry coefficients a 3ym (A) and 
fsym(A) shown as a function of mass number for three tem- 
peratures. The full and dashed lines refer to parametrization 
I and II. The black, blue and red lines correspond to T —0, 
4, and 8 MeV, respectively. 
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FIG. 3: (color online) The thermal dependence of the coef- 
ficients /3e and /3f (from parametrization I) is displayed in 
panels (a) and (b) . The thermal evolution of the surface sym- 
metry coefficients a s aym and fs ym (from parametrization II) is 
shown in the right panels (c) and (d). The blue and red lines 
in both the left and right panels refer to calculations with the 
KDEO and SBM interactions, respectively. 



T = 0, the value of j3 E or /3 F is 12.1 and 13.9 MeV for the 
SBM and KDEO interactions, respectively; they compare 
well with the value of ~ 13 MeV obtained from analyses 
of the 'experimental' symmetry energies of isobaric nuclei 
@- With temperature, /3e decreases for both the inter- 
actions; @f shows a nominal increase. We have, however, 
noticed that both a v sym / '/3e and fs ym /PF are nearly tem- 
perature independent, lying in the range of ~ 2.64±0.01. 

The temperature-dependent surface symmetry coeffi- 
cients a s sym and f% ym as used in Eqs. (2) and (11) in 
the definition II are shown in the right panels of Fig. 3. 
At T = 0, a s sym is 44.8 MeV and 39.2 MeV for the 
KDEO and SBM interactions, respectively, close to the 
phenomenological value of ~ 45 MeV [||. With temper- 
ature, for both the interactions, a s Bym increases sharply 
showing the growing importance of the surface term in 
a sym {A). The surface free energy coefficient fl ym , how- 
ever, displays a slow decrease with temperature for the 
KDEO interaction. As for the SBM interaction, fg ym is 
nearly temperature-independent. 

A comparison with calculations in Rcf. [l|| may now 
be in order. In both calculations, the surface symme- 
try coefficient seems to be more sensitive to temperature 
compared to the volume symmetry coefficient. However, 
in Ref. [l6| . in the limited temperature range they ex- 
plore, the temperature dependence of the surface coeffi- 
cients seem to be more pronounced than those seen in 
the present calculation. There are subtle differences too, 
the lack of self-consistency of the density profiles used 

• FT^l _l _ 1 i i l_* l_ J 
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FIG. 5: (color online) The thermal dependence of the symme- 
try coefficients of nuclei shown for three mass number. The 
full and dashed lines refer to parametrization I and II. The 
black, blue and red lines correspond to A —60, 140 and 220, 
respectively. 

proximations involved may be the reason behind these 
differences. 

In Fig. 4, the mass dependence of the a sym (A) and 
fsym{A) is shown at three temperatures, T = 0, 4 and 
8 MeV. Panels (a) and (b) in the figure display a sym (A) 
for KDE0 and SBM interactions, respectively; panels (c) 
and (d) display f sym (A). The full lines correspond to 
definition I for the symmetry coefficients, the dashed lines 
do the same for definition II. The general findings are : 
for a particular mass number, a sym (A) decreases with 
temperature, f S ym{A) increases. At fixed temperature, 
asym(A) and f sym {A) increase with A; this follows from 
the definitions. The values of f sym {A) seem to depend 
little on the parametrization I or II; similar is the case 
for a sym (A) except at very high temperature. 

For the fixed values of nuclear masses, the tempera- 
ture dependence of a sym (A) and f sym (A) are exhibited in 



Fig. 5. The masses chosen are A = 60, 140 and 220. Pan- 
els (a) and (b) in this figure display a sym (A) for KDE0 
and SBM interactions, respectively; panels (c) and (d) 
do the same for f sym (A) . The black lines pertain to A = 
60, the blue lines to A = 140 and the red lines corre- 
spond to A = 220. The general findings in Fig. 4 that 
&sym{A) falls and f S ym{A) shows a very slow increase 
with temperature is reinforced from this figure. For the 
SBM interaction, a near constancy of f sym (A) with a 
slight dip in the middle of the temperature range is seen. 
This was also occasionally observed earlier [23( with a 
different definition of f sym (A) - in the spirit of Eq. (3). 
It is also observed that for a chosen interaction, both 
parametrization I and II yield nearly the same value of 
the symmetry coefficients except for A = 60 at higher 
temperatures. 

To summarize, in a liquid-drop-modcl-inspired fit of 
the total energies and free energies of a system of nuclei 
evaluated in a subtraction-implemented finite tempera- 
ture Thomas-Fermi framework, the temperature depen- 
dence of the symmetry energy coefficients of nuclei have 
been evaluated in this communication. Two different en- 
ergy density functionals, one with the zero-range Skyrme- 
type KDE0 and the other with a finite-range SBM inter- 
action have been employed for this purpose. The general 
behavior of the temperature dependence of the symmetry 
coefficients seems to be nearly independent of the energy 
functional used. For cold systems, the calculated vol- 
ume and surface symmetry energy coefficients lie within 
the constraints set from analyses of different experimen- 
tal data. With temperature, the symmetry free energy 
coefficients show a weak change. A strong temperature 
dependence of a v sym is however observed, the tempera- 
ture dependence of a s sym is even stronger; this results in 
a rapid fall in a sym {A) of the atomic nucleus as the tem- 
perature rises. The calculations, in addition throw light 
on the thermal mapping of the volume and surface en- 
ergies which are in excellent qualitative agreement with 
those in common usage. 
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